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Abstract 

I construct the spectral function of the Luther-Emery model which describes 
one-dimensional fermions with one gapless and one gapped degree of freedom, 
i.e. superconductors and Peierls and Mott insulators, by using symmetries, 
relations to other models, and known limits. Depending on the relative mag- 
nitudes of the charge and spin velocities, and on whether a charge or a spin 
gap is present, I find spectral functions differing in the number of singulari- 
ties and presence or absence of anomalous dimensions of fermion operators. 
I find, for a Peierls system, one singularity with anomalous dimension and 
one finite maximum; for a superconductor two singularities with anomalous 
dimensions; and for a Mott insulator one or two singularities without anoma- 
lous dimension. In addition, there are strong shadow bands. I generalize the 
construction to arbitrary dynamical multi-particle correlation functions. The 
main aspects of this work are in agreement with numerical and Bethe Ansatz 
calculations by others. I also discuss the application to photoemission exper- 
iments on ID Mott insulators and on the normal state of ID Peierls systems, 
and propose the Luther-Emery model as the generic description of ID charge 
density wave systems with important electronic correlations. 

PACS numbers: 71.27.+a, 71.30.+h, 71.45.Lr, 79.60.-i 
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I. MOTIVATION 



Non-Fermi liquid behavior in correlated fermion systems is an exciting topic of current 
research. One-dimensional (ID) correlated electrons (more precisely: one-dimensional quan- 
tum systems with gapless excitations) are a paradigmatic example of non- Fermi liquids: their 
low-energy excitations are not quasi-particles but rather collective charge and spin density 
fluctuations which obey each to their proper dynamicsB. The key features of these "Lut- 
tinger liquids"! are (i) anomalous dimensions of operators producing correlation functions 
with non-universal power-laws, parametrized by one renormalized coupling constant K v per 
degree of freedom v = p (charge), a (spin) which have the status of the Landau parameters 
familiar from Fermi liquid theory; (ii) charge-spin separation, leading to a fractionization 
of an electron into charged, spinless, and neutral, spin-carrying collective excitations, with 
different dynamics determined by velocities v p ^ v a . Each of these features leads to (iii) ab- 
sence of fermionic quasi-particles. Responsible are the electron-electron interaction which is 
marginal in one dimension and therefore transfers nonvanishing momentum in scattering pro- 
cesses at all energy scales, and the nesting properties of the ID Fermi surface. They produce 
divergent 2k p charge and spin density fluctuations which then interfere with Cooper-type 
superconducting fluctuations. 

All three features clearly show up in the single-particle spectral functionlhi 

p(q, uj) = — 7r -1 Im G(kp + q, p + uj) (1.1) 

which can be measured (within the "sudden approximation" ) by angle-resolved photoemis- 
sion (ARPES) [with bad angular resolution, one essentially measures N(u) = J2 q p(<L w ) an d 
is able to probe only features (i) and (iii)]. The spectral function is purely incoherent! - !, 
at best with peaks at the dispersion energies of the elementary charge and spin excitations, 
indicating that the electron behaves as a composite particle built on more elementary exci- 
tations. In Eq. ( |1 . lp , G is the Fourier transform of the retarded electronic Green's function 

G(xt) = -i6(t)({*(arf),* t (00)}) , (1.2) 

kp the Fermi wave number, and p (= 0) is the chemical potential. 

Much experimental effort has been devoted to studying and attempting to "prove" Lut- 
tinger liquid correlations in various quasi- ID systems. Examples are organic conductors of 
the family based on the molecule TMTSF (Bechgaard salts) where both NMPii and (par- 
tially) photoemissionB have provided evidence in favor of a Luttinger liquid picture, quantum 
wires fabricated into semiconductor nanostructuresl, or edge states in the fractional quan- 
tum Hall effect!. In all cases however, there appear to be problems with the precise values of 
the parameter K p derived, or with some other aspects of the interpretation in terms of a Lut- 
tinger liquid. It is not clear to date to what extent these discrepancies are due to the neglect 
of some experimentally important factor in the theory (such as, e.g. three-dimensionality or 
electron-phonon coupling in the chain systems, or deviations from the special filling factors 
in the quantum Hall edge states), or indicative of more fundamental problems either with 
theory or experiment. 

ID (organic and inorganic) charge density wave (CDW) systems apparently could provide 
an alternative field of search for these typically one- dimensional correlations. Photoemission 
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indeed has produced results^ similar to the Bechgaard salts when performed with low angu- 
lar resolution. With high angular resolution, a broad dispersing feature has been identified 
in (TaS 64)2 M while two such signals have been measured in the blue bronze K . 3 MoO 3 K 
Even though the actual situation in K Q3 Mo0 3 may be slightly more complictated because 
there are two almost degenerate bands cutting the Fermi energy, it is clearly of impor- 
tance to first understand the photoemission spectrum expected from the metallic phase of a 
single-band CDW material. Finally, while this paper was prepared, new experiments on the 
organic two-chain conductor TTF — TCNQ became available which clearly show dispersing 
signals both on the TTF and TCNQ chains with very unusual lineshapesS. Specifically, 
the TCNQ signals are somewhat similar to K 03 MoO 3 , and we know from independent ex- 
periments that there are strong 2&^-CDW fluctuations on this chain in the metallic stateB. 
(The TTF-chain exhibits strong Akp-CDW fluctuations at very high temperature and is 
expected to be a Luttinger liquid.) 

The association of the two dispersing signals of K 03 MoO 3 with the charge and spin 
excitations of a Luttinger liquid is suggestive. As I will explain in the next section in more 
detail, it is incompatible, however, with the CDW transitions observed in these materials. 
This incompatibility motivates the consideration of the Luther- Emery model and is at the 
origin of the work reported here. Section II will discuss this model, its generic role as a 
low-energy fixed-point of ID quantum systems which have both gapped and gapless degrees 
of freedom, and the picture we had of its correlations prior to this work. 

Recently, photoemission experiments also have been performed on the ID Mott insulator 
SrCuoM- In Mott insulators, the charge fluctuations are gapped while the spins remain 
gapless. Their low-energy physics, therefore, can again be described by a Luther-Emery 
model, and our theory can be adapted to study the spectral functions of ID Mott insulators. 
Earlier, angle-integrated photoemission on BaVS 3 has been interpreted as evidence for a 
Luttinger liquidS. The behavior of the conductivity, however, is more insulator-like, and 
present theory might be of interest there, too. 

Section III presents the construction of the single-particle spectral function ( |1.1| ). In 
Section IV I present results for the spectral functions of the spin-gapped Luther-Emery 
model, i.e. ID Peierls systems and superconductors. In Section [V], the spectral functions 
of ID Mott insulators are presented. Section [Vj shows how the construction procedure 
of Section [III] can be generalized to arbitrary correlation functions of local operators. I 
compare my results with information from other studies in Section [VI 1] and use them for an 
interpretation of published experiments in Section |V111| . I conclude with a short summary 
and a brief perspective. Partial results have been presented earlierEH 



II. THE LUTHER-EMERY MODEL 

The Luther-Emery model extends the Luttinger model by including the backscattering 
of electrons across the Fermi surface. Its Hamiltonian isEl 

H = H + H l + H 2 + H i , (2.1) 
H o=J2 v F(rk - k F ) : cl ks c rks : , (2.2) 

r,k,s 
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H i = E [dn 5 *,*' + 9i±Ss,-s'] fdx : : ( 2 - 3 ) 

#2 = y E [^2||^, S ' + #2iA,- s '] p+,s{p)p-, s '{-p) , (2.4) 

^4 = E [^4||^,s' + fi , 4lA,-s' : Pr,sip)Pr,s'(-p) ■ ■ (2.5) 

c r fc s describes fermions with momentum and spin s on the two branches (r = ±) of the 
dispersion varying linearly [e r (k) = vp{rk — about the two Fermi points ±kp, ^ ryS (x) 
is its Fourier transform, and 

PrM = E : C t,k+p,s C r,k,s : = E (4,fc+p,,Cr,fc, S - ^,0<4 I fc,« c r,fc,»)o) ( 2 -6) 
fe fe 

is the density fluctuation operator which obeys a bosonic algebra 

[Pr,s{p),Pr's'{-p')] = -8r,r'8s,a>Sp,p>-7r- ■ (2-7) 

The Luttinger model is obtained for gi = and includes only forward scattering. 

In one dimension, fermions can be transformed into bosons, and for the Luttinger model, 
there is an exact operator identity relating a fermion operator ty rs (x) to the bosonic density 



ct or> 

fluctuations (|2.(j|)Eil. For our purposes, the approximate expression 



ty rs (x) ~ lim -== exp -= [r$ p (x) - G p (x) + s {r$ a (x) - Q a (x)}\ . (2.8) 
y2na \y2 ) 

with the two phase fields 

?7T „-a\p\/2-ipx 

^)-yE- [u + (p)+u.(p)] , (2.9) 

and 

j-jr „-a\p\/2-ipx 

e,(^yE- [i/ + (p) - i/_(p)] , (2.10) 

found earlier by Luther and Peschefi is sufficient. 

This formula allows for a boson representation of the Hamiltonian and of all correlation 
functions. Before, it is important, however, to recall the physics of the phase fields <& v (x) 
and Q u (x) in ( [2.8| )EEI1. The charge density fluctuation operator is related to & p (x) by 
J2 r Pr( x ) = _ n~ 1 d^ p (x)/dx, and likewise for spin a. When an additional particle is inserted 
into the system, a kink of amplitude n is formed in $> u (x). These fields therefore describe 
the scattering phase shifts of the particles present in the system, generated by the particles 
added. The operators inserting the particles are exponentials of the dual fields Q v (x) = 
J U v (x)dx, where Yi v is the momentum conjugate to $„: [IL^x), $ v {x')\ = —i5(x — x'). In a 
general fluctuation operator whose correlation function we wish to evaluate, the prefactor of 
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iQ u / a/2 measures the number of //-particles it inserts into the system while the prefactor of 
i<& u / \/2 measures the number of //-particles it rearranges at constant total //-particle number 
to generate the desired fluctuation. By //-particle, we label, in the first place, 



m ru {x) = (27ra)- 1/4 exp{-i[r<5> u (x) - B v {x)]/y/2} , 



(2.11) 



the slowly- varying charge or spin part of the fermion operators ty rx (x) but, with phase 
factors reflecting the appropriate Fermi seas, these particles will describe the holons and the 
spinons of the ID Bethe-Ansatz soluble models. 

The boson form of the Luther-Emery Hamiltonian becomes 



Hq + H 4 = — ( nV F + 9iu) ■ U r {p)u r (-p) 

urpj^O 

+#2 = 7 XX 2 ^ - 9i\\)v+(p)v-(-p) , 

L up 



H 



2g 



dx cos 



X 



v r (p) are the operators for the charge and spin densities 

Pr(p) = ~J= [PrM +Pr,lip)] 



(2.12) 
(2.13) 
(2.14) 



and the interactions have been transformed as 



9i 



(2.15) 



(2.16) 



Diagonalizing the Luttinger part (i.e. H excluding Hi±) generates the renormalized veloci- 
ties of the collective charge and spin excitations and their stiffness constants 



\ 



Vf + 



9au 

7T 



92u - gin/2 



TV 



7TVF + g4u ~ 92u + glJ|/2 
F + 94u + 92u -01||/2 



(2.17) 



The phase fields transform as Q v {x) — > $ v (x)^K u and Q u (x) — > Q v {x) j \[K^. The main 
effect of the (74-interaction is a renormalization of v v . We therefore drop if 4 from explicit 
consideration in the following, and always assume correctly renormalized velocities v v . 

For K a — 1 sufficiently large with respect to backscattering is irrelevant, and the 

Luther-Emery model reduces to a Luttinger liquid. Its renormalized value of K a can be 
calculated, e.g. by perturbative renormalization groupil which is well-controlled in this 
case or, if applicable, fixed to unity by the requirement of spin-rotation invariance. Charge 
and spin excitations are gapless, and depending on the value of K p , the dominant are spin 
density wave (SDW, K p < 1, repulsive forward scattering) or triplet pairing (TS, K p > 1, 
attractive forward scattering). Charge density wave (CDW) and singlet superconducting 
(SS) fluctuations, respectively, are sub dominant. 
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The backscattering Hamiltonian H%x is, for K a — 1 small enough compared to \g\s\ 
relevant perturbation and opens a gap A CT in the spin excitation spectrum 



e a (q) = ±y/vtf + Al . (2.18) 

Luther and Emery have shown that for the special value K a = 1/2, the interaction Hamilto- 
nian Hi± ( p.!4|) can be represented as a bilinear in spinless fermions, using the bosonization 
formula (|2.8j ) for spinless fermions (multiply the argument of the exponential by y/2 and 
drop the ex-fields), and diagonalizedEj. On this Luther- Emery line K a = 1/2, the gap is 



computed exactly to be A CT = \gi±\/27ra [a is an infinitesimal in ( |2.8j ) but often associated 
with a cutoff of the order of a lattice constant] . Renormalization group then allows to derive 
the gap for arbitrary K a . The charges remain gapless. 

The Mott insulator is the consequence of an instability in the charge channel, caused 
by Umklapp scattering off the lattice for commensurate band-fillings. The Umklapp Hamil- 
tonian appropriate for a half-filled band is obtained by simply replacing spin by charge in 
Eq. ( p.!4| ), and its coupling constant often is denoted by g s ±. Here the spins are gapless 
while relevant Umklapp scattering opens a gap A p in the charge channel. This generic pic- 
ture applies (with little modification only) to all even commensurabilities (k F a = [r/s]n/2, 
s even). The situation is different for s odd, where the Umklapp operator necessarily cou- 
ples charges and spinsS, and we exclude these cases from our study. The Mott insulator is 
dominated by 4/c^-CDW and/or SDW correlations. 

While the Luther- Emery solution is essentiallyii exact, it is useless for computing cor- 
relation functions since there is no practical relation between the physical fermions and 
the spinless pseudofermions. Still, we have some qualitative information on the correlation 
functions. Several methods" support the idea that, in the gapped phase, correlations of 
the $ CT -field tend towards a non-zero constant as \x\ or \t\ — > oo while those involving ex- 
ponentials of its dual field Q a (x) decay exponentially in space (or oscillate in time). The 
spin gap quenches low-energy spin fluctuations, therefore SDW and TS correlations should 
be exponentially suppressed. With a constant asymptotic value of § a , CDW and SS are 
enhanced with respect to a Luttinger liquid, and now dominate over SDW and TS. The 
opening of a spin gap is a necessary condition for the emergence of dominant SS or CDW 
correlations in a ID metal. As a corollary, a Luther-Emery phase must exist in the normal 
state of CDW systems (or superconductors) between a Luttinger liquid and the 3D ordered 
low-temperatures phases. One therefore should be careful in interpreting the properties of 
the metallic "normal state" of a CDW system (or of a ID superconductor) in terms of a 
Luttinger liquid. 

For the one- and two-particle spectral functions, there is a general belief that the opening 
of a gap affects the system for frequencies smaller than this gap while the behavior of the 
ungapped system is essentially recovered at larger frequency scales. The exponential decay 
(resp. oscillations) of correlation functions involving operators exp[z(. . would cut off 
(shift) the divergences as functions of q (u) they had possessed in the Luttinger model. 
Possibly important power-law prefactors to exponentials have not been discussed. There 
has been almost no calculation or systematic construction of such functions - in particular 
dynamical onesil - and, to mvknowledge, no critical check of these hypotheses by numerical 
simulation prior to this worklHM. 
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A wide variety of models fall into the Luther-Emery universality class and my results 
should be applicable there in a low-energy sector: Luttinger liquids coupled to phonons and 
related models so long as they are incommensurate, have wide regions of parameter space 
with gapped spin fluctuations and gapless charge the negative-?/ Hubbard model at any 
band-filling has a spin gapS, and the positive-f/ Hubbard model at half-filling has a charge 
gapHil; with longer-range interactions, charge gaps can occur at different rational band- 
fillings, too. The t — J-model has a spin gap at low density!!. Spin gaps occur freauently in 
models of two Luttinger or Hubbard chains coupled by single-particle tunnelingHS. Also 
when a 2fcp-CDW is established in many coupled Luttinger chains as a consequence of inter- 
chain Coulomb interaction, the system passes through a region of attractive backscattering 
which opens a spin gapil. 



III. CONSTRUCTION OF THE SPECTRAL FUNCTION 

I now present a systematic construction of the single-particle spectral function, Eqs. 
( |1.1D and ( |1.2j ), for the spin-gapped Luther- Emery model. The Green's function exhibits the 
full complexity of the problem, involving all four phase fields G^, while many others 
are easierS. They will be discussed in Section |VT[ Here, we limit ourselves to the diagonal 
terms of the Green's function, both in the branch index r and in the spin index s, and 
further assume spin-rotation invariance, so that s is dropped alltogether. This assumption, 
which I will make throughout the paper unless exceptions for the sake of an argument are 
stated explicitly, further implies K a = 1. With the nonvanishing expectation values of 
operators exp[z(. . .)$o-] generated by the gap opening, finite off-diagonal terms are possible, 
in principle, both here and in multi-particle correlation functions. They can be calculated 
in complete analogy to the terms discussed here, and we ignore them in the following. 

Using bosonization (|2.S| ), the retarded Green's function for right-moving fermions (r = +) 
can be represented product 

G(xt) = -iQ(t)e lkFX [g p (xt)g a (xt) + (x -»• -x, -t)} , (3.1) 

of charge and spin correlation functions 

g v {xt) = (V +U (xt)*l u (00)) . (3.2) 

The product structure is a consequence of the charge-spin separation of the Hamiltonian 



T]). The spectral function (|1.1|) then is a convolution 



P(<1, u) = ttA? / dq dJ [g p (q',uj')g a (q - q',u - uj') + (q -> -q , u -> -uj)} . (3.3) 



(27T) 



The charge part is easy and can be calculated in the Luttinger model (I only display the 
leading uj- and g-dependence) 

g p (q, u)~e(uj- v p q)Q(uj + v p q)(uj - v p q)^\u + v p q)^ 1 ' 2 (3.4) 
with 7 P = {K p + Kp 1 — 2)/8 for K p ^l , 



Q(u + v p q) 
^uT+lj~q 



6{uj - v p q) {K p = 1) . (3.5) 
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Using a similar expression for the spins, one can reproduce in detail the spectral func- 
tions of the Luttinger model calculated elsewhere directlyHl. Notice that the divergences 
are stronger than for a spinless Luttinger model ensuring that singularities remain after 
performing the convolution integrals. For both K v ^ 1, the coalescence of three of the 
four singularities of g p (q,u) and g a (q,u) is needed to generate a singularity in the spectral 
function of the Luttinger model; if one of them, e.g. K a , is unity, the coalescence of two 
singularities is sufficient. 

The determination of the spin correlation function is more involved because it has no 
simple representation in terms of the Luther-Emery pseudofermions, excluding any exact 
calculation. I now show that the leading behavior of this function can, however, be uniquely 
constructed from symmetries, equivalences, and known limits if the Ansatz is made that 
g a {xt) is a product of power laws and exponentials in x and t. There is a variety of arguments 
requiring this form, and we will give them in the following, together with the construction 
procedure. 

The important steps are: (i) Representing the Hamiltonian in terms of right- and left- 
moving fermions requires g a to be a function of x ± v a t only. In general, g a will contain both 
power laws (/±) and exponentials (/exp) of these variables 

9a{xt) ~ f + (x - v a t)f_(x + v a t)f cxp (x 2 - v\t 2 ) . (3.6) 

Interactions other than g 4 can only mix left- and right-moving excitations, producing prod- 
ucts of x ± v a t, or functions thereof, but cannot introduce new dependences on x and/or t. 
This is consistent both with the boson solution of the massless Luttinger phase and with 
the Luther-Emery solution of the gapped phase. (The Lorentz invariance of the Luther- 
Emery model requires all correlation functions of Luther-Emery pseudofermions to depend 
on x 2 — v 2 t 2 only - and by implication all those of the physical fermions whose operators 
can be represented in terms on Luther- Emery fermions alone.) The exponential part f exp 
necessarily is a function of x 2 — v 2 t 2 only. All dependences on x and t other than through 
functions of x 2 — v 2 t 2 must therefore be present also in the Luttinger model (gi± = 0), and 
necessarily are of power-law form. 

(ii) The limit of a vanishing gap A CT — > can also be used to constrain the function g a (xt), 
but is rather subtle. To make the argument clear, we momentarily relax the assumption 
of spin rotation invariance so that the spin channel of the model is described by gi± and 
general K a . (Alternatively, we can look at a Mott problem with Umklapp scattering gs± 
and K p ^ 1 is more natural.) In the limit A a — > 0, the function f exp (x 2 — v 2 t 2 ) — > 1 here, 
because the exponential dependences are introduced by the finite gap. Straightforwardly, 
one would now identify the product f + {x — v a t)f^(x + v a t) = CT (xt)\I/J.(00)) Slx =o with the 
spin part of the spectral function of the remaining Luttinger model, i.e. Eq. ( |3.7| ) below 
with anomalous exponents <5_ = [K a + K^ 1 — 2)/8 and 5 + = 5- + 1/2. This physically 
appealing procedure was used in an earlier paperEl, and possibly could describe the physics 
of a small-gap Luther-Emery model. 

Taking the limit A CT — > to constrain eventual power-laws in g a (xt) involves different 
physics, however, and the above argument must be modified. For vanishing gap, g a must 
reduce to the correlation function of the free Luttinger model (K a = 1), no matter what 
value of K a would describe the hypothetical Luttinger model obtained from the Luther- 
Emery model ( |2.1| ) for gu_ = 0, i.e. independently of any assumption on spin-rotation 
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invariance. Physically, this is so because the anomalous operator dimensions K a ^ 1 of the 
Luttinger model are a consequence of singular low-energy virtual particle-hole excitations. 
When there is a gap at the Fermi surface, these processes are quenched, and one is left with 
the exponent K a = 1 of the free mo defi Notice that this argument implies that we consider 
a rather large gap. 

Accidentally, the spectral functions given earliei0 remain correct. This, however, is due 
to the limitation to spin-rotation invariant interactions there. They impose K a = 1 for the 
power-law functions f±(x =F v a t) in any case. 

With f exp (x 2 — vlt 2 ;A a = 0) ~ 1 one can determine all possible power-laws f± up to 
corrections varying more slowly than a power law, to be 

f±(x =F v a t) = [a - i(x =F v a t)]~ s± (3.7) 

with exponents 

5+ = 1/2 , 6- = . (3.8) 

These are the exponents of a free Luttinger correlation function for the spin part of a 
right-moving fermion. I remphasize that they arise because of the quenching of low-energy 
particle-hole excitations by the spin gap and hold independent of any assumption on spin- 
rotation invariance. (As we will see below, the corresponding result for the charge channel 
implies that there cannot be any anomalous dimensions in a ID Mott insulator with spin- 
rotation invariance respected). 

(hi) From the equivalence of the Luther-Emery model to a classical 2D Coulomb gasil 
(using the Matsubara formalism of imaginary times r = it, putting y = v a r) and Debye 
screening of the charges above the Kosterlitz-Thouless temperature, one deduces an expo- 
nential factor 



f cxp (x ± v t) ~ exp(-cA a ^/x 2 /v 2 ~t 2 ) (3.9) 

with an undetermined constant c, in / exp . This equivalence quite generally excludes any 
decay faster than 



In this picture, the perturbation Hamiltonian ( 2.14 ) generates a Coulomb gas of charges 
q e = ±1, and the ^-fields of the Green's function appear as two test charges q' e = ±1/2 
whose (bare logarithmic) interaction is modified by screening from the Coulomb gas. The 
gapped Luther-Emery phase corresponds to the high-temperature plasma phase of un- 
bound charges in the Coulomb gas, and the screening can then be treated in the Debye- 
Hiickel approximationil. Here, the effective potential between the charges is V(r) ~ 
exp(—K£)\r\)/ \J n D \r\ with the Debye wavevector k d = 2A CT /t> cr il. The CT -fields can then 
be viewed as magnetic monopoles with strengths q m = ±1/2. Their interaction is again 
logarithmic, and they couple to the electric charges with V em (f) ~ — arctan(y/x)il. Clearly, 
the high-temperature plasma of electric charges e e = ±1 modifies the effective monopole- 
monopole interaction which becomes 

^=-^ + Tf-) 2 -JTT ' ( 3 - 10 ) 
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where I have used the Debye-Hiickel polarization propagator 



m = fS^ • (3-ii) 

Fourier-transforming back to real space, one obtains 

Kn- m (r) ~ In \f\ + c'\f\K D (3.12) 

with an open constant d oc c. One observes an antiscreening effect here: in the presence 
of the electric charges, the magnetic monopoles are confined more strongly than without 
charges! Going back to real times, ( |3.12| ) produces the exponential dependence in ( |3.9| ) and, 
most importantly, gives additional (in fact, for those multi-particle correlation functions 
which only depend on x 2 — v 2 t 2 the only firm) justification for the presence of power-law 
prefactors in addition to exponential terms in (|3.6|). 

(iv) The open constant c in (|3^ ) can be determined from a spectral representation of / exp , 
and our interpretation of the bosonization formula ( [2.8|) . Fourier transforming f exp (x, t), one 
obtains 

Q(lu 2 - v 2 a 2 - c 2 A 2 ) 
= 2^cA (T ( J_^f_ c2A , ) -; 2 (3.13) 

which has a gap of magnitude cA a in its spectrum. This gap must correspond to the excita- 
tion of \n\ spinons where n is the prefactor of iQ a (x)/\/ r 2 in the operator whose correlation 
function we wish to calculate. This constrains the prefactor in the exponential to c = |n| 
quite generally. For the single-particle Green's function n = 1, and we obtain c = 1 here. 

(v) The present construction of g a (xt) is not an exact calculation. It is therefore im- 
portant to look for exactly known cases which can be used as tests, to confirm the validity 
of this construction. Gulacsi has calculated explicitly the t = O-Green's function of a ID 
Mott insulator!!: He finds G(x) ~ exp(— A p |x|)/|x| which is in complete agreement with 
the present theory when the 1 / ^/jxj-contribution from the ungapped channel is multiplied 
to Eq. ( 3.14|) below. That there may be a power-law prefactor in the charge part of the 



spectral function has also been realized but well hidden in publications, by others^. 

In Section 0, I will discuss further tests of these rules based on two-particle correlation 
functions. 

From the rules (i) - (v), I find 

g a (x,t) ~ exp ^-A^yx 2 - tffi/vA /\/a + i(v a t - x) . (3.14) 
Fourier transformation then gives 



The comparison of (|3.15 ) with ( |3.5| ) (after p — > a there) is interesting. The 5-function 



translates the absence of anomalous dimensions in the gapped channel, a consequence of rule 
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(ii), rather than spin-rotation invariance as in the a- version of ( |3.5|) . The change in dispersion 
due to the spin gap enters through this 5-function. The frequency-dependent prefactor is 
the same as in the gapless system. However, due to the different argument in the 5-function, 
it no longer becomes singular in the limit q, uj — > but has an upper limit of A" 1 / 2 now. 
A similar effect occurs in the Green's function of ID quantum antiferromagnets, where the 
opening of the spin gap cuts off a singularity of the prefactor of the delta functionEH The 
factor in parentheses is a coherence factor translating the enhanced spin-pairing tendency at 
the origin of the spin gap, and one readliy recognizes the same structure as for the coherence 
factors u q , v q familiar from the theory of superconductivity. 



IV. SPECTRAL FUNCTION FOR THE SPIN-GAPPED LUTHER-EMERY 

MODEL 



We now must convolute g a (q, uj), Eq. ( |3.15| ), with the charge part, Eqs. ( |3.4j ) or (|3.5|). The 



results depend on the relative magnitudes of the charge and spin velocities. We therefore 
treat separately the cases of (A) repulsive interactions (in the sense that the effective forward 
scattering matrix element g2 P — <7i||/2 > 0), i.e. K p < 1, and v p > v a , where Peierls-type 
2/c^-CDW fluctuations dominate, and (B) attractive forward scattering, i.e. K p > 1 and 
v p < v a , when singlet superconducting fluctuations are most important. (The inequalities 
on the velocities and K p usually go with each other as listed when standard lattice models 
are treated. Of course, when one takes all as free parameters, other combinations are 
possible. Relevant for the subsequent classification then are the velocities.) 

What could we expect from our knowledge of the Luttinger liquidi? There the singulari- 
ties at uj = v p f a )q arise from processes where the charge (spin) contributes all of the electron's 
momentum q and the spin (charge) none. The same argument applied to the Luther-Emery 
model predicts signals at the renormalized spin dispersion e a (q), Eq. Q2.1S| ), and at a shifted 
charge dispersion 

e P (q) =v p q + A a . (4.1) 

Figure 1 shows the location of the signals expected from this argument. The A^-shift in 
the charge dispersion comes from the fact that the zero-momentum spin fluctuation can 
only be excited at a cost of A CT . As will be seen below, however, the spectral functions of 
the Luther-Emery model never show two singularities with these dispersions. The intuitive 
predictions on the spectral function of the Luther-Emery model basically transcribe the 
standard argument that the behavior of correlation functions is modified on energy scales 
below the gap (correlations are suppressed there) but recovered almost unchanged on higher 
energy scales. Our results will show that for dynamical, q and cu-dependent correlations, 
this argument is not trustworthy. 



A. One-dimensional Peierls "insulators" 



We assume v p > v a and K p < 1, implying dominant CDW correlations. Calling these 
systems "insulators" is a misnomer, however, because the charges are gapless and the systems 
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are metallic. More precisely, we think about the Luther-Emery model here as describing the 
"normal" metallic state above a CDW transition. 

The convolution of g a and g p , Eq. ( |3.3| ), is rather straightforward now. After executing 
the a/-integral, singularities are obtained from the coalescence of the two singularities carried 
by g p (q,u). The result of the calculation is shown schematically in Fig. 2 for q < (unlike 
previous papers, we present the spectral functions as those of the occupied states, i.e. as 
they would be measured by a photoemission experiment). There are indeed features at the 
special frequencies shown in Fig. 1. On the spin dispersion e a (q), there is a true singularity 

p[q,u w -sM)] ~ - e*b)][-«> ~ ^M)T~ 1/2 (4-2) 

as in the Luttinger model. Here, a is defined as a = (K p + K^ 1 — 2)/4 = 2j p since the notion 
of a K a does not make sense in a spin-gapped system. Folklore would then predict another 
singularity \to + ^((jOI^ 1-1 ^ 2 (short dashed lines in Fig. 2) which is not observed here. It is 
cut off instead to a finite maximum of order 

p[q,un-e p {q)]~^r m ■ (4-3) 

The reason for cutting of the Luttinger divergence on the charge dispersion is related to the 
non-singular prefactor (for q — > 0) in g a (q,u), cf. Eq. ( [3.15|) and the subsequent discussion, 



and the convolution makes this effect apparent on the charge dispersion e p (q). The spin 
gap therefore supresses the divergence associated with the charge dispersion while on the 
renormalized spin dispersion, the spectral response remains singular. 

At positive frequencies, the Luther-Emery model has pronounced shadow bands. Here, 
the Luttinger liquid only has very small weight. The weight in the Luther- Emery model is 
much stronger, and the spectral function has the same overall shape as at negative frequen- 
cies. For q < 0, the negative frequency part is enhanced by a coherence factor 1 — v a q/e a (q) 
while a factor 1 + v a q/e a (q) decreases its shadow. These factors translate the increased 
coherence due to the spin pairing and the finite spin gap, and are a consequence of the 
corresponding coherence factors in Eq. Q3.15| ). Of course, as suggested by Fig. 1, one can 
also view the shadow bands as bending back from the Fermi (or more precisely: the gap) 
energy as k is increased beyond hp. This view perhaps is closer to a real photoemission 
experiment. 



B. One-dimensional superconductors 

We now take v p < v a , i.e. attractive forward scattering. This implies K p > 1, and such a 
system has dominant singlet pairing fluctuations. Interestingly, two true singularities occur 
here whose location is shown in Fig. 3. There is one singularity on the renormalized spin 
dispersion 

p[q,u « -e a {q)\ ~ ©(? " <7c)e[-w - e ff (q)][-u - £ CT (g)r 1/2 

+ e(q c -q)\- UJ -e a (q)r 1 / 2 , (4.4) 

which is one-sided for q > q c and two-sided for q < q c . 
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is a critical wave vector which arises in the convolution procedure from searching the mini- 
mum of s a {q') + v p (q — q') as a function of q' . At this wavevector, the dispersion 

e P (q) =e a (q c ) +v p (q-q c ) (4.6) 

is tangential to e a (q). For q < q c , a divergence 

p[q,uj ^ -e p (q)] ^ Q(-uj - e p (q))[-uj - i p (q)}^ (4.7) 

on this shifted charge dispersion splits off the spin divergence. Again, there are strong 
shadow bands with the same functional forms as the main bands, specifically with two 
singularities, and with intensities controlled by coherence factors. The dispersions of the 
signals are displayed in Figure 3, and the shape of the spectral function is sketched in 
Figure 4. 

Notice that, quite generally, that the behavior of p(q, to ~ ±A CT ) is determined by that of 
the spin part close to A CT and that of the charge part at cu « 0. Unlike earlier conjectures^, 
it is therefore not necessary to know details of the charge dynamics on a scale to ~ A a where 
the Luttinger description may have acquired significant corrrections. 

The fc-integrated density of states then is N(u) ~ Q(u — lA^^o; — |A CT |) a , independent 
of the magnitudes of the velocities. There is no weight below the gap, and the typical gap 
singularity in the density of states of the spin fluctuations is wiped out by convoluting with 
the gapless charges. 

It is quite clear now that certain properties of ID fermions - the dynamical ones involving 
(1+1)D Fourier transforms - are affected by the gap opening on all energy scales, contrary 
to common expectation, while those depending on one variable alone are modified only on 
scales below the gap energy. Despite the opening of a gap in the spin channel, singular 
spectral response remains possible in q- and a»-dependent correlation functions. 



V. SPECTRAL FUNCTION OF ONE-DIMENSIONAL MOTT INSULATORS 

The spectral function of a ID Mott insulator can be computed as a special case of the 
generic solution presented above. One simply has to change a «-> p everywhere and put 
K a = 1 in the gapless spin channel for spin-rotation invariance (which we assume to hold, 
again). Importantly, the exchange of p and a also applies to the inequalities on the velocities 
v u , where again two cases must be distinguished. 

Both factors g u in the convolution now involve 5-functions. In the case of repulsive 
forward scattering v p > v a , one now finds a spectral function with two singularities, similar 
to the case of a ID superconductor. Since K a = 1, the anomalous single-particle exponent 
a = 0, i.e. one obtains two inverse square-root singularities. In the main band (u < for 
q < 0), the spectral function becomes 

/ \ Q(g-gc)Q[-^-e P (g)] , Q(q c - q)Q[-uj - e a (g)] 
p(q,uj) ~ = 1 — (5.1) 

y/-u-e p (q) yJ[-LJ - e a (q)]\ -u-e p (q)\ 
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with £ a (q) = s p (q c ) + v a (q — q c ). An important difference to the case of a superconductor 
occurs in the shadow band: since the spectral function of the gapless spin channel has 
no shadow band of its own, the singularity on e a (q) in the shadow band is missing. The 
shadow band therefore has a single singularity on the charge dispersion e p (q) with a weight 
depressed by a coherence factor with respect to the weight of the main band signals. The 
effect is completely analogous to the appearence of a single nonanalyticity in the (very weak) 
shadow bands of a Luttinger liquid with spin-rotation invariant interactionsiHl. The shape 
of this spectral function is sketched in Figure 5. The location of the singularities follows 
Figure 3 with the replacement p <-> o except for the shadow bands where the straight lines 
should be ignored. 

The case v a > v p again is different. Compared to the case of the ID Peierls "insulator", 
the anomalous dimension a on the charge dispersion drops out due to spin-rotation invari- 
ance, giving an inverse square-root singularity on e p (q). Also the finite maximum on the 
shifted spin- dispersion e a {q) does not occur. This is because the <5-function has zero weight 
in the energy domain where the square-root prefactor in Eq. (|3.15|) takes its maximum. The 
shadow band, of course, has a single inverse-square-root singularity with the usual coherence 
factors. Thus, the spectral function for this case becomes 

K^)~ e ' M ," £ f' ■ (5.2) 

up to coherence factors, and the density of states 

N{uS) ~ Q(\u\ — A p ) x regular function . (5.3) 

The spectral properties of a doped Mott insulator, of course, depend on the detailed sce- 
nario emerging from a more complete theory. Work on the Hubbard model shows, however, 
that the upper Hubbard band qualitatively survives a finite dopant concentration!!!^ Con- 
tinuity then suggests that as the insulating state is left by varying the band-filling, spectral 
weight is gradually taken out of both the main and shadow bands of a spectral function 
such as those discussed before, and transferred into the charge and spin divergences of a 
Luttinger liquid signal. Although the spins are left unaffected in the transition and only 
a charge gap opens, both the charge and the spin signals are predicted to be shifted and 
strongly modified by doping. This is a direct consequence of the convolution property ( |3 .3| ) 
of the single-particle spectral function. When superposing (to a first approximation) the 
two signals, care must be taken, in addition, to account for the dependence of the chemical 
potential on doping level. 



VI. GENERALIZATION TO OTHER CORRELATION FUNCTIONS 

We now discuss the construction of other correlation functions for the Luther-Emery 
model. Clearly, due to charge-spin separation, they can again be written as convolutions of 
charge and spin correlation functions. Consider a general local operator 

0^\x)=^{x)^ v {x) , (6.1) 
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where ^ rv {x) had been introduced in Eq. (|2.11|) , and a positive (negative) exponent is 



understood as a creation (annihilation) operator. Bosonizing V1 the "Iv-field acquires a 
prefactor (m—n), and Q u is multiplied by (m+n) with respect to the single-particle operator 
^ru- If gapless channel is assumed to be the charge v = p (as we have done throughout this 
paper except in the preceding section), the correlation function of O p behaves as 

R^ n \xt) = (o^ n \xt) [O( m ' n) (00) m 

,2; 



/ Luttinger 

/ i*i • \—in 2 /2 i | • i ■ • \— n 2 /2 

{a + tVpt — ix) 1 (a + iv pt + ix) ' 



(m — n) 2 / T r 1 * (m+n) 2 / T r— 1 

x U+*v) 2 +*v ' (6 ' 2) 

Its Fourier transform is 

R<-™'"Xq,u) ~ 0(u> - K p9 )e(a> + t, p ?)(a> - ty?) 7 '"""'"'^ + v^y'"'^ , (6.3) 
lf-> = f + ^=^Vp - 1) + " 1 

7 " " 2 + 8 8 

We now turn to such an operator for spin, CT , in the presence of a spin gap. When 
the spin gap opens due to the Hamiltonian ( |2.14| ), the <3> CT -field develops long-range order. 
Its dual field, G CT , then is disordered, and its correlations will contain exponential terms 
similar to / exp , Eq. Q3.9|). We now have to distinguish two cases, (i) If we have m = —n, the 
operator 0(. m '~ m ) can be represented in terms of the ^-field alone. Since this is the ordering 
field, we simply can put it to a constant value, implying R^ 1 '" 1 ^ (xt) ~ 1, and the space- 
time dependence of the total correlation function is then determined by the charge part 
R^^^q^uj) (which may carry different powers s,t of ^/ r ,p, depending on the spin directions) 
alone, and given by Eq. ( |6.3| ). One can, in principle, go one step further and account for 
the long- wavelength fluctuations out of the ground state- value of Q a . A convenient method 
for this again is the mapping onto a classical 2D Coulomb gas. Since the $ CT -field of the 
correlation functions introduces electric test charges, we know that in the massive Luther- 
Emery phase their interaction is exponentially screened (cf. Section p|) . We then find the 
fluctuation contribution 



9 exp (-2Kdx*/vl-tA 
<[$^)-M00)] 2 > \ x y v ^_ t i )1/ , J • (M 

I will discuss an interesting application in a moment. 

However, if (ii) m ^ —n, the spin correlations contain the disorder field 0^ dual to the 
$0-, and the gap opening will lead to exponential factors as in Eq. ( |3.9| ). This is the case 
for the Green's function, cf. Eq. ( |3.2|) . We apply the same rules (i) - (v) as in Section 
ll|. Specifically, the prefactor of the gap in the exponential is c = \m + n|, by comparing 
the energy for the insertion of |m + n\ a-particles into the system with the gap obtained in 
the spectral representation of the exponential. The power-law prefactor is that of the free 
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Luttinger model because there cannot be any anomalous dimensions in a gapped fermion 
system. In (xt)-space, the correlation function then is 

R^ n \xt) = (o^ n \xt) [O( m '")(00) lf 



~ (a + iv a t — ix) m 7/2 (a + iv a t + ix) n ^ 2 exp ^— \m + n\/\ a \J x 2 / v 2 — t 2 ^ . (6.5) 

This expression can be Fourier transformed and convoluted with an appropriate charge part. 

What the present construction cannot do, however, is to give information on the mag- 
nitude, or a possible vanishing, of the prefactor of the correlation function. One example is 
the 2fc^-CDW correlation function in the half-filled replusive Hubbard model, where a naive 
use of the construction above would predict (in real space at t — 0) a dependence ~ x _1 
which, on physical grounds, is not expected to be important in that modefi Qualitative 
information can be obtained in that situation from renormalization group studies, where 
one can monitor how the amplitude of a correlation function changes as one moves away 
from a Luttinger liquid fixed points. A complete solution of this problem presumably would 
require an exact boson representation of the physical fermions in a Luther-Emery model, 
including fermion raising operators. 

To conclude this Section, I discuss two more test cases for my construction procedure. 
Consider the transverse 2/cp-spin-correlation functions©'!! 

R S Dw±(xt) = (O S Dw±(xt)Ol DW± (00)) = (* t _ i (^)^ +T (^)^ t +>T (00)*_ i (00)) (6.6) 

in the Luther-Emery spin-gap regime. The spin density wave operator can also be repre- 
sented as 

S DW±(xt) = ^^expj-zv^ [$„(x) + Q a (x)]} = e^O^ix^O^^xt) . (6.7) 

We now limit ourselves to the spin component of the correlation function and obtain, using 
Eq. (pi), 



exp {-2& a Jx 2 /vl-t*) 

R^- l \xt) V , V . (6.8) 

'x 2 — v 2 t 2 



Fourier transformation gives 



jfr^-rrtl- 4 *?- («■») 



Iuj 2 - v 2 a q 2 - 4A^ 

One the other hand, on the Luther-Emery line K a = 1/2, one can refermionize the operator 

0^- x) {x) = V2wVt(x)tf + (x) (6.10) 
in terms of spinless fermions by inverting the spinless variant^ of the bosonization 



formula ( |2.8| ). The limitation of this procedure to the Luther- Emery line is inessential 
because different coupling constants will only affect the magnitude of the spin gap but not 
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the form of the excitation spectrum, so long as A CT > 0. Now, one can calculate Ry 1 '" 1 ' (q, lo) 
as the pairing correlation function of spinless fermions in a fermion representation. Such a 
calculation has been outlined by LeeE3, and the result derived from his expressions agrees 
with Eq. ( |6.9| ) both concerning the regions of nonvanishing spectral weight, and the critical 
exponents of the singularities. Incidentally, my own expressions are more complicated than 
Lee's by additional terms and additional occupation functions n(k) and 1 — n(k). They 
conspire with the coherence factors [l±v a q/e a (q)} to produce a prefactor v\q 2 / (t>^g 2 +4A^) to 
the leading inverse-square-root singularity which vanishes asg^O. At g = 0, a subleading 
term oc 0(|u;| — 2A CT ) times a regular function remains. Apart these subtle prefactors, the 
exact fermionic calculation reproduces the result of the construction procedure advocated 
here for the correlation functions of the Luther-Emery model. 

A final test is provided by the charge correlations of a ID Mott insulator. In general, 
the charge density operator n(x) has contributions at wavevectors q ~ 0, 2kp, 4kF, etc. 



y/2d$ p (x 

n(x) ~ 



7r dx ixa 
2 

(2tto) 



+ — exp \-2ik F x + \/2i®Jx)} cos[v / 2$ (T (x)] 



+ exp [-Aik F x + V8i$ p (x)} . (6.11) 



In a half-filled band, = 2ir/a, a reciprocal lattice vector so that the 4fcp-term effectively 
does not oscillate when measured on the lattice sites. When the Mott gap A p opens, the field 
$ p orders at a finite constant value. The third term in ( p. lip then translates the long-range 



charge order, the first term measures the long-wavelength fluctuations out of this ordered 
ground state, and the second term measures 2kp charge fluctuations. Using the arguments 
at the beginning of this section (after a <->• p), we obtain from the first two terms a spectral 
function 



q MU 2 2 V } **} ] • (6.12) 



The zero- frequency 5-function comes from the "4fcir"-part, and the high-frequency signal 
from the d$ p /dx-term. In principle, one could also calculate the 2fci?-part. However, ex- 
perience with the Hubbard model suggests that prefactors not specified here suppress the 
2fcp-CDW fluctuations on the lattice sitesS, and we do not consider them here (similar, 
and nonvanishing contributions, however appear in 2fcp-SDW correlation functions, and in 
a "bond order wave" which is best described as a 2/cj?-CDW centered midway between two 
sites). 

The spectral function R n (q,u) has been calculated recently by Mori and FukuyaniE 
They do not give an explicit expression which would allow to check the critical exponents, but 
the region of nonvanishing spectral weight, and the overall shape of the high-frequency signal 
are consistent with Eq. ( 6.12Q , whereas the 5-function in Eq. ( 6.12j) seems to be missing. It 



is present, however, in a numerical diagonalization of an extended Hubbard model^, and 
provides another, though more superficial test of our construction. 
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VII. RELATION TO OTHER WORK 



In the preceding sections, we have discussed some tests for the dynamical correlation 
functions of the Luther-Emery model constructed hereHS Independent verification comes 
from work on many models which fall into the Luther-Emery universality class. 

In particular, numerical studies have attempted to look into the spectral properties of cor- 
related fermion models. Quantum Monte Carlo simulations of the ID Hubbard model at half- 
filling, a prototypical Mott insulator with v p > v a provides evidence for pronounced shadow 
bands, much stronger than those of the doped systems which form Luttinger liquidsB At 
present, the resolution is not good enough to directly visualize the two dispersing inverse- 
square-root singularities found here. However, recent improvements on doped Hubbard 
mo delS lend hope that Quantum Monte Carlo will be able, in the near future, to confirm 
the predicitons made here. 

The ID t — J-model at half-filling also forms a Mott insulator with v p > v a , and exact 
diagonalization of lattices up to 22 sites has allowed a calculation of the spectral function of 
this modelEl. While the location of regions of finite spectral weight, and of the singularities 
agrees with the present study, numerical diagonalization on such small systems does not 
allow to determine the critical exponents of the divergences of the ID Mott insulator. 

Spin gaps also arise in many lattice models. E.g. for two coupled Luttinger, Hubbard, 
or t — J-chains, there are wide regions of parameter space where the spin fluctuations are 
massive, and the single-particle spectral function has been calculated occasionally^. Again, 
exact diagonalization finds important shadow bandsi but the resolution is not good enough 
to separate the two dispersing divergences found in Section [TV| for a superconductor, not to 
speak of the much weaker signal on the shifted charge dispersion e p (q) predicted above for 
a CDW system. 

Evidence for such a weak signal, and for a divergent signal on the gapped spin dispersion 
e a {q) comes, however, from exact diagonalization of a t — J — J'-model where a spin gap 
opens for certain values of J'e3. These authors observe a very strong spinon signal, the holon 
peak is anomalously weak, as predicted here. 

A Bethe Ansatz calculation of spectral functions for a ID Mott insulator has recently 
been performed by Sorella and Parola (SP) based on the ID supersymmetric t — J modeled, 
and also confirms essential aspects of the present work. In their model, v p < v a so that 
we predict a single inverse-square-root singularity on e p (q). Such a singularity is also found 
from the Bethe Ansatz solution used by SP. When a finite magnetization is included, SP find 
critical exponents which explicitly depend on the momentum of the hole created. One would 
expect from universality and the possibility to transform a positive-f/ Hubbard model into 
one with negative U by a particle-hole transformation on one spin species alone, that such 
spectral functions should also describe spin-gapped systems with v p > v a . We do not find 
such momentum-dependences in the work presented here. SP's method, however, requires 
the calculation of the ground state and low-energy properties of the spin Hamiltonian at 
finite total momentum of the spin system. These explicitly depend on the momentum, and 
produce the momentum-dependent exponents. In the Luther-Emery model, one calculates 
a spinon excitation with some momentum with respect to a zero momentum ground state. 
The momentum-dependent correlation exponents found by SP certainly are beyond scope 
and possibilities of the present model. On the other hand, their method does not allow to 
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look into more subtle features than critical exponents, such as the finite maximum which 
we found in this case. 



VIII. APPLICATIONS TO EXPERIMENTS 

Importantly, our results could prove useful in the description of the photoemission prop- 
erties of certain quasi-lD materials. 

There have been angle-resolved photoemission experiments on the ID Mott insulator 
StCu02 with a gap 2Ap ~ 1.8e The lineshapes observed were anomalously broad and 
showed unsual dispersion. As a consequence, the authors proposed a description in terms of 
a system with charge-spin separation, where the broad feature would, in fact, be composed 
of the unresolved spin and charge signals. In addition, a strong shadow band bends back 
from the gap edge for k > kp. Its dispersion is consistent with the one of the charge signal 
for k < kp. Clearly, these observations are fully consistent with the theory presented here, 
which predicts two inverse-square-root singularities beyond some critical wave vector (cf. 
Fig. 5), and a single one below, as are the accompanying diagonalization results on a ID 
t - J-modeM 

More interesting in the present context are a number of unexplained ARPES results 
on organic and inorganic materials which undergo Peierls transitions at low temperatures. 
Specifically, ARPES experiments on the blue bronze K0.3M0O3 by several groups show two 
dispersing peaksB Also in the organic conductor TTF — TCNQ, anomalous lineshapes are 
observedli3. Of interest here is the TCNQ-band which shows 2/c^-CDW fluctuations in the 
metallic state@ and triggers a series of transitions into a low-temperature CDW phase. While 
some materials such as the Bechgaard salts!, or the TTF-band of TTF — TCNQ (which has 
strong 4fcp-CDW fluctuations") may well fall into the Luttinger liquid universality class, 
it is particularly surprising that CDW systems such as K , 3 MoOs, or the TC AQ-band in 
TTF — TCNQ, should behave as Luttinger liquids. In fact, the photoemission properties 
are in striking contrast to the established picture of a fluctuating Peierls insulator which has 
been applied quite universally to describe the normal state of CDW systems^. It predicts a 
strongly temperature dependent, narrow [|u;| < A C dw(T = 0)] pseudogap and p(q < O^jj) 



is governed by a broadened quasi-particle peak at uo < and a weak shadow at uo > 

A Luttinger liquid interpretation for the CDW photoemission is highly suggestive but 
encounters problems which are all resolved in a Luther-Emery framework, (i) As has been 
explained before, Luttinger liquids have no dominant 2A;^-CDW correlations: for repulsive 
interactions (K p < 1), spin density waves are logarithmically stronger than CDWs§, and the 
behavior of lattice models is consistent with this pictured. For attractive interactions, the 
system is dominated by superconductivity^. A spin gap is a necessary condition for promot- 
ing CDW correlations in correlated ID electron systems and is realized in the Luther- Emery 
model! (ii) 2fcp-CDWs often are due to electron-phonon coupling, and renormalization group 
provides us with a detailed scenari The dependence of the spin gap on electron-phonon 
coupling A, the phonon frequency wd, and K p , can be calculated reliablyS A spin gap 
also opens if 2/c^-CDWs are caused by Coulomb interaction between chains!!, (iii) The 
spin susceptibility of CDW systems above the Peierls temperature decreases significantly 
with decreasing temperature indicative of activated spin fluctuations. This applies applies 
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both to K0.3M0O3 at temperatures from Tp to beyond 700 K0, and to the TCiVQ-chain 
in TTF — TCNQ where the magnetic susceptibility contributions of both chains can be 
separated by NMPiIll. Notice in this context that at finite temperature, the density of states 
in the spin channel of the Luther-Emery model is essentially the same as for the Lee-Rice- 
Anderson theory of a fluctuating Peierls insulatoril, implying that both models will have 
similar x{T). The temperature-dependent susceptibility alone therefore cannot discriminate 
between these two theories. Remarkably however, in Kq^MoO^ the conductivity is metal- 
lic in the same temperature range: early experiments over a restricted temperature range 
find the resistance p(T) ~ tH while very recent data taken to much higher temperatures 
even suggest a sublinear temperature dependence!! - not unlike the one found in Luttinger 
liquids with repulsive electron-electron interactions!!. In TTF — TCNQ, p{T) ~ T has 
been foundil, but it is not known how the individual chains contribute to this dependence. 
The experiments are incompatible with the temperature dependence of the conductivity 
expected in a fluctuating Peierls insulatorEl which indeed is observed in some organic ma- 
terials and also {TaSe^I ■ (iv) For a Luttinger model, the stronger divergence in p(q,u) 
is associated with the charge mode and disperses more quickly than the weaker signal. In 
the experiment on K0.3M0O3, the quickly dispersing signal is less peaked than the slow one. 
On the other hand, the important feature of the Luther-Emery spectral function, Fig. 1, is 
that the spin gap supresses the divergence of the charge signal which disperses more quickly 
than the divergent spin contribution, (v) A CDW transition out of a Luther-Emery liquid 
by opening a charge gap at the Peierls temperature, is also consistent with subtle transfers 
of spectral weight in regions away from the Fermi energy, observed in spectra taken through 
the true CDW transition^. In these experiments, the spectral weight at the Fermi energy 
is essentially zero at any temperature. However, at some finite energy below Ep, the weight 
drops with a temperature dependence consistent with a BCS-like gap. In a naive charge-spin 
separating, Luther-Emery scenario, one would postulate the opening of a charge gap A p at 
the Peierls temperature (as a consequence of the establishment of 3D coherence, allowing 
for the finite-T transition), in addition to the preexisting spin gap. Thus one expects a 
drop of spectral weight at the Peierls transition in an energy range between Ep — A a and 
E F — A a — A p which, on a sufficiently coarse temperature scale, would amount to a shift 
of the leading edge by A p . More likely, the establishment of 3D coherence will destroy to 
some extent the ideal spin-charge separation of the ID Luther-Emery model, and produce 
a single CDW gap Aqdw > A a below the transition, both for charges and spins. 

On a quantitative level, there is one major problem for the description of the normal state 
of most CDW systems: the spin gap derived from an analysis of the magnetic susceptibilities 
is much smaller than the spin gaps derived from the peak maxima of the ARPES signals. At 
present it is not clear if this indicates a fundamental problem with a Luther-Emery model 
(the problem would however not be solved with any competing theory), if this is due to some 
ununderstood effect in the photoemission process, or if it is due to some extrinsic sample 
property. In another language, it is not clear what mechanism is responsible for apparent 
gaps which systematically are a sizable fraction of the valence band widths. 

This phenomenology is not consistent with many other theories proposed for ID fermions. 
Theories based on a fluctuating Peierls insulator would have to explain the two dispersing 
bands seen in Kq. 3M0O3 as two separate bands. Two such bands indeed exist but the 
implication would be that band structure calculations get one of them too narrow by a 
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factor of 5, but get the correct dispersion for the other oneS. Moreover, they cannot 
reconcile the activated susceptibility with the essentially metallic conductivity above the 
Peierls temperature. 

Standard Luttinger liquidffl, but also the anomalous ground states obtained from cou- 
pling Luttinger chains so long as their low-energy fixed point is a Fermi liquidil, do neither 
produce the CDW correlations, nor the activated susceptibility Notice, however, that both 
transversely coupled Luttinger liquids (Kopietz et a/.S ) and the ID Hubbard model@ can, 
under some circumstances, produce spectral functions where the peak on the spinon dis- 
persion is stronger than that on the charge dispersion. They, however, would predict a 
Fermi surface crossing of the photoemission signal which is not observed experimentally, in 
addition to the problems listed above. In the experiments, instead, the dispersing spectral 
features bend back from the Fermi energy as k is increased beyond kp, in a manner strongly 
reminiscent of the shadow bands discussed before. 

Depsite (important) quantitative problems, the Luther-Emery spectral function is con- 
sistent with the photoemission experiments on Kq^MoO^ and TTF — TCNQ, and beyond 
that, the model is consistent with much of the other experimental phenomenology available. 
I emphasize that while the agreement of the Luther-Emery spectral function with the ob- 
served photoemission lineshapes certainly is an argument in favor of this model, it is the 
consistency of its predictions with most other experiments available which suggests that it 
might be a natural starting point for a description of the low-energy physics of these CDW 
materials. 

Obviously, this suggestion is somewhat speculative and independent support is called 
for. Its virtue is that it comes to grips with the puzzle that the spin susceptibilities of 
K Q 3 Mo03 and TTF — TCNQ decrease with decreasing temperature while the conductivity 
are metallic, that it leaves space for the good description of optical properties as a fluctuating 
Peierls insulator (they only probe the charge fluctuations which will form CDW precursors 
at temperatures much below the spin gap opening, presumably as a consequence of emerg- 
ing 3D coherence), and that it provides an (admittedly phenomenological) description of 
the photoemission properties of these materials with extremely ID electronic properties!!. 
As in the Bechgaard salts!, a single-particle exponent a ~ 1/2... 1 would be required 
implying strong long-range electron-electron interactions, and there is at best preliminary 
support from transport measurement sS, for such strong correlations in Kn^MoO^. Retarded 
electron-phonon coupling could increase a over its purely electronic valueEl. To what extent 
this mechanism contributes could be gauged from the measured a which must be larger than 
the one derived from the enhancement of v p over the band velocity (helas strongly depending 
on the accuracy of band structure calculations). In TTF — TCNQ, the analysis is made 
difficult by the presence of two chains. There is evidence for strong long-range electron- 
electron interactions on the TTF-chain from the observation of Akp-CDW fluctuations, but 
the situation for TCNQ is less clear. If a sizable enhancement of the dispersion of the 
ARPES signals over the estimated bandwidths can be interpreted as evidence for long-range 
electronic correlations, they would indeed be present on both chains. 
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IX. SUMMARY 



In this paper, I have presented a construction of the dynamical correlation functions 
of the ID Luther-Emery model. This model has one gapped degree of freedom, and an 
ungapped one, and describes ID superconductors and Peierls insulators (spin gap) and ID 
Mott insulators (charge gap). It is a natural extension of Luttinger liquid theory to the 
peculiar phase intermediate between metal and band insulator, made possible in one di- 
mension by the phenomenon of charge-spin separation. The dynamical correlation functions 
presented here show where and to what extent the two typically ID features of a Luttinger 
liquid: charge-spin separation, and anomalous dimensions of operators, survive in the pres- 
ence of a gap in one channel. Since an exact calculation of such correlation functions usually 
is not possible in a Luther-Emery model, our construction relied heavily on limiting cases, 
symmetries, and equivalences to other models. However, it successfully passed several tests 
in situations where exact results were available from other methods. 

The main emphasis of the paper was on the single-particle spectral function which is mea- 
sured in photoemission. We showed that, generically, charge-spin separation and anomalous 
dimensions are also visible in the spectral functions of the Luther- Emery model. Specifically, 
for a spin gapped system with repulsive interactions, describing a ID charge density wave 
system, the spectral function has a true singularity on the gapped spin dispersion with an 
anomalous exponent a — 1/2 while on the charge dispersion, the Luttinger liquid divergence 
is cut off to a finite maximum by the spin gap - a results which finds a straightforward 
explanation in terms of convolution of charge and spin correlation functions. For attrac- 
tive interactions, i.e. a ID superconductor, two divergences with anomalous dimensions are 
found. For ID Mott insulators, i.e. a charge gap, one finds one or two inverse-square-root 
singularities, i.e. no anomalous dimension (due to spin-rotation invariance), depending on 
the order of the velocities of the charge or spin fluctuations. It was also shown how these 
procedures can be generalized to two- and multi-particle correlation functions. 

Besides predicting spectral functions for the many ID models falling into the Luther- 
Emery universality class, there are a few experimental situations where these results can 
be usefully applied. They successfully describe the photoemission spectrum of the ID Mott 
insulator SrCuO^, to an extent leaving few questions open, the most notable one being 
experimental resolution. Less clearcut but perhaps more interesting are CDW materials such 
as K0.3M0O3 and TTF — TCNQ which show very unusual photoemission spectra. These 
are qualitatively consistent with a Luther-Emery model, and we have proposed that these 
materials might, most naturally, be described in this framework. A Luther-Emery phase is 
necessary as an intermediate between a Luttinger liquid and a long-range ordered CDW, and 
K 3 MoO3 and TTF — TCNQ are natural candidates for searching for such a strange metal. 
This scenario requires strong electron-electron interactions at least at high energies, and not 
all CDW materials need fall into this scheme. If the electron-phonon interaction is so strong 
as to produce CDW precursor fluctuations at very high temperature, and the electronic 
correlations are weak enough, the establishment of a Luttinger liquid, and the crossover to 
a Luther-Emery liquid at lower temperature, may be quenched, and a fluctuating Peierls 
insulator!! or a bipolaron liquid^ maybe a more appropriate picture. Some CDW materials 
such as (TaSe 4 ) 2 M, (perylene) 2 PF 6 lll, (fluoranthene) 2 PF(@ apparently are consistent with 
this picture. However, K 3 MoO 3 and TTF — TCNQ are not consistent, and the consistency 
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of the spectral functions constructed in this paper with the published experiments, and the 
analysis of further experiments indicate that, besides electron-phonon coupling, electronic 
correlations must be important in these CDW systems. 
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FIGURES 



FIG. 1 Dispersion of peaks in the spectral function p(q, oj) of a spin-gapped Luther- 
Emery model with v p > v a . The dispersion laws e p (q) and e a (q) are given in the text. The 
heavy solid and dashed lines give the signals in the main band [sign(u;) = sign(g)] while the 
light dashed lines label the shadow bands [sign(u;) = - sign(g)]. 

FIG. 2 Spectral function of the spin-gapped Luther-Emery model for q < 0. v p > v a has 
been assumed, as applies to a ID Peierls insulator. The dashed line at —e p (k) indicates the 
Luttinger liquid divergence which is suppressed here to a finite maximum. 

FIG. 3 Dispersion of singularities in the spectral function p(q, uj) of a spin-gapped Luther- 
Emery model with v p < v a . The solid lines give the signals in the main band while the dashed 
lines label the shadow bands. 

FIG. 4 Spectral function p(q, uj) of the spin-gapped Luther-Emery model with v p < v a 
(applying to ID superconductors) for q < 0. 

FIG. 5 Spectral function p(q, uj) of the charge-gapped Luther- Emery model, describing 
ID Mott insulators, with v p > v a for q < 0. The dispersions of the signals follow Fig. 3 with 
p <-> a everywhere, and the straight dashed lines in the shadow bands must be ignored. 
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